Measuring the decoherence rate in a semiconductor charge qubit 
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We describe a method by which the decoherence time of a solid state qubit may be measured. 
The qubit is coded in the orbital degree of freedom of a single electron bound to a pair of donor 
impurities in a semiconductor host. The qubit is manipulated by adiabatically varying an external 
electric field. We show that, by measuring the total probability of a successful qubit rotation as 
a function of the control field parameters, the decoherence rate may be determined. We estimate 
various system parameters, including the decoherence rates due to electromagnetic fluctuations and 
acoustic phonons. We flnd that, for reasonable physical parameters, the experiment is possible with 
existing technology. In particular, the use of adiabatic control fields implies that the experiment 
can be performed with control electronics with a time resolution of tens of nanoseconds. 

PACS numbers: 03.65.Yz, 03.67.Lx, 73. 21. La, 63.20.Kr 



I. INTRODUCTION. 

It is generally believed that scalable quantum comput- 
ing devices will eventually be fabricated in solid state 
systems, and many ideas have been proposediSiSi^i^ 
Of these, only superconducting systems have, as yet, 
yielded single qubit devices capable of demonstrating a 
large number of single qubit rotationsi^ No multi qubit 
devices have been reported in the solid state. A single 
qubit device however is extremely useful as it enables an 
experimental measurement of the qubit decoherence time 
to be made. This number will ultimately determine if a 
particular solid state implementation is scalable (that is, 
capable of reaching the error threshold required for fault 
tolerant operation^). Decoherence refers to the uncon- 
trollable coupling between the degree of freedom coding 
the qubit and other degrees of freedom in the qubit 's en- 
vironment. Such uncontrollable interactions lead to the 
qubit becoming entangled with these inaccessible degrees 
of freedom, with the result that the state of the qubit is 
not precisely defined by its preparation and subsequent 
control by unitary gates. Under such circumstances, the 
outcomes of direct measurements on the qubit are de- 
scribed by a mixed state, corresponding to an average 
over the inaccessible degrees of freedom. 

In solid state systems the sources of decoherence are 
legion, and include phonons, nuclear spins, and electro- 
magnetic fluctuations. Which sources of decoherence are 
relevant depend on what particular degrees of freedom 
are used to encode the qubit. A great deal of experi- 
mental and theoretical work remains to be done if we are 
to achieve understanding of the limitations of solid state 
implementations of qubits. In this paper we will focus 
on one particular qubit encoding based on the electron 
charge degree of freedom. (Decoherence of the charge de- 
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FIG. 1: A schematic representation of the double donor, single 
electron system. The qubit is encoded as the spatial localiza- 
tion of the electron charge, relative to the donor sites. The 
surface gates may be used to control the bias term, e, in the 
qubit Hamiltonian. The radio frequency single electron tran- 
sistor (RF-SET) may be used to readout the position of the 
electron. 



gree of freedom is also of relevance to schemes in which 
spin qubits are coupled by an exchange interaction^*2i& 
since charge decoherence can lead to leakage errors dur- 
ing exchange interaction gates>i£) 

To be specific, we will consider a system which consists 
of two phosphorus donors, embedded in a silicon sub- 
strate, which share a single excess electroniii The device 
is depicted schematically in Fig. The qubit is coded 
in terms of the relative position of the electron. We de- 
note localized single particle states by {\L), \R)} where 
\L) corresponds to an electron localized on the left donor 
site, while \R) denotes an electron localized on the right 
donor. These states are not eigenstates of the Hamil- 
tonian when the potential is perfectly symmetrical. We 
may however represent localized states in terms of the 
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two lowest states of the potential; the symmetric ground 
state \Eg) and the antisymmetric first excited state \Eas), 
by \L,R) ^ {\Es) ± \Eas))/V2. The system prepared in 
one of these localized states will oscillate coherently be- 
tween them at the tunnelling frequency A = Eas — Eg. If 
the potential is biased sufficiently far from symmetry (by 
applying an external electric field) the localized states 
become good approximations to the energy eigenstates. 

The single particle Hamiltonian for the double donor 
system may be approximated by the two level qubit 
Hamiltonian 

^^qb = -h^^CFz - fiyCTj: , (1) 

where cr^ = \L){R\ + \R){L\ and = \L){L\ - \R){R\ 
are the Pauli spin operators in the basis, and 

e{t) describes the bias of the potential away from sym- 
metry, due to an external electric field. We have indi- 
cated that this may be a function of time, as discussed 
below. For a nonzero bias, the energy gap between 
the instanta neous grou nd state and first excited state is 
hE{t) = hy/A^+e{ty. 

This approximate Hamiltonian is valid if A, |e|, 
ksT/h <C wqi, where loqi is the angular frequency cor- 
responding to transitions between the ground and first 
excited states of an electron bound to a single isolated 
donor. For phosphorus donors in silicon, wqi = 1.8 x 10^'^ 
rad s~^»i^ As we discuss in Sec. |nl the tunnelling fre- 
quency, A, depends on the distance between the donors. 
For a donor separation of around 40 nm, the tunnelling 
frequency is approximately A ^ 10^" rad s^^. 

Decoherence in this system can be due to phonons that 
cause transitions between the energy eigenstates of the 
system. As we show in Appendix 1X1 however, the cor- 
responding timescale for such transitions can be made 
much longer than all other timescales in the problem 
by choosing an appropriate donor separation. Interac- 
tions with electromagnetic fluctuations in the environ- 
ment (e.g. due to thermal voltage noise in nearby sur- 
face gates) however is more serious. In this paper we will 
model such processes using the spin boson model. This 
model has been extensively discussed in the literature 
(see for example Refs. fisl and fl^ . 

Our objective is to find a way to experimentally de- 
termine the decoherence rate. It might be thought that 
this is easily done by monitoring the decay of the co- 
herent tunnelling oscillation, by allowing the system to 
evolve for a time t and then determining the expected po- 
sition of the electron relative to the double donor system, 
{az){t). Repeating for a number of different values of t, 
and observing the decay time of the oscillations in (a^) (t) 
would yield the decoherence rate. While this is possible 
in principle it is difficult in practice, because the coher- 
ent evolution must be turned on and off (for example, by 
rapidly changing the bias field s{t)), on timescales much 
shorter than the reciprocal of the tunnelling frequency, 
E{t)~^ . Using this technique, the tunnelling frequency 
itself must be much larger than the decoherence time. 



which is expected to be of the order of nanoseconds (see 
Sec. II V Bll . Therefore, measuring the decay of coherent 
oscillations directly would require accurate switching of 
the qubit Hamiltonian on a timescale of tens of picosec- 
onds. Despite these difficulties, a similar experiment has 
been achieved in a superconducting charge qubit i 

Other experiments have focused on continuous mea- 
surement of the charge degree of freedom of excess elec- 
trons in a coupled quantum dot system, using a nearby 
quantum point contact electrometer A signature of 
charge decoherence in the coupled dot system was ob- 
served by monitoring the average current through the 
electrometer, although a large contribution to the ob- 
served decoherence rate is thought to be due to the back 
action of the electrometer on the coupled dot system. 
This back action is due to the shot noise of the electrons 
tunnelling through the quantum point contact. 

In a recent paper^ii an alternative method was pro- 
posed to determine the decoherence rate for flux qubits 
implemented in a radio frequency SQUID systemii^ 
Rather than attempting to observe the decay of coher- 
ent oscillations of the flux, the authors proposed that the 
qubit polarization be reversed by adiabatically sweeping 
the qubit Hamiltonian parameters. They argued that 
the decoherence time can be determined by observing the 
probability of success of the adiabatic inversion process 
as a function of the parameter sweep time. 

In this paper, we describe a scheme for determining 
the decoherence rate in the single electron, double donor 
system described above. Our scheme also makes use of 
adiabatic manipulation of the Hamiltonian parameters. 
We show that an experimental estimate of the decoher- 
ence rate can be obtained by preparing the system in the 
ground state under strong positive bias (a state locafized 
on the left donor), adiabatically sweeping the bias to zero 
{e{t) — 0) and then holding the bias at zero for a pe- 
riod thoid, before adiabatically sweeping to the opposite 
bias and then determining whether or not the system has 
changed its localized charge state. The final charge state 
of the system can be measured using a radio frequency 
single electron transistor (RF-SET),i2iSa An RF-SET can 
be kept in a quiescent state during the qubit evolution, 
and therefore the detector back action should not add a 
significant contribution to the observed decoherence rate. 
A plot of the probability of finding the electron on the 
right donor site versus ihoU will in general fall from a 
value close to unity, to substantially less than unity, over 
a timescale determined by the decoherence rate. 

The advantages of this method over one in which coher- 
ent oscillations are directly observed are twofold. Firstly, 
substantially fewer measurements are required, since it 
is not necessary to plot out several coherent oscillations. 
Secondly, the timescales over which e{t) must be var- 
ied are determined by the decoherence timescale itself, 
rather than the (much shorter) timescale for coherent os- 
cillations, E{t)~^ . Estimates of the relevant parameters, 
presented below, suggest that the experiment can be per- 
formed with electronics with a time resolution of tens of 
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FIG. 2: Approximate energy gap, A, between the lowest sym- 
metric and antisymmetric eigenstates of double donor, single 
electron system, under zero bias, as a function of donor sep- 
aration. 



nanoseconds, rather than tens of picoseconds. 

This paper is organized as follows. In Sec. |n]we esti- 
mate the tunnelling frequency, A, for the double donor, 
single electron system, as a function of the donor separa- 
tion. In Sec. mil we describe the scheme for determining 
the decoherence rate in more detail. In Sec. IIV Al intro- 
duce the spin-boson model for the coupling of the qubit to 
the environment. In Sec. IIV Bl we calculate an estimate 
for the strength of the system-environment coupling for 
the case of decoherence due to thermal voltage noise in 
nearby surface gates, and present the results of numeri- 
cal calculations of the evolution of the qubit under such 
a coupling. In order for the experiment to be viable, 
a number of constraints must be satisfied. We quanti- 
tatively discuss these in Sec. |3 and also find a set of 
experimentally achievable parameters that satisfy these 
constraints. We also discuss a number of other issues 
relating to the implementation of this scheme. 



II. APPROXIMATE ENERGY LEVELS OF THE 
SINGLE ELECTRON, DOUBLE DONOR SYSTEM 

The tunnelling frequency. A, may be estimated by de- 
termining approximate energy eigenvalues for the low- 
est energy symmetric and antisymmetric eigenstates for 
the double donor, single electron system. Finding ex- 
act values for these energy levels is complicated by the 
fact that the conduction band electron dispersion rela- 
tion in silicon is anisotropic, and also by the valley-orbit 
interactionr— For the purposes of this work, however, 
it will be sufficient to gain an order of magnitude esti- 



mate for A. To this end, we ignore the conduction band 
anisotropy and assume that the localized states \L) and 
|_R) may be represented by Is-orbitals centered on the 
left and right donor sites respectively. We take the Bohr 
radius for an isolated donor state to be 

flfj = esi ( — ) as , (2) 

where esi = 11.7 is the dielectric constant for silicon, 
TOg is the mass of a free electron, rriT = 0.2me is the 
transverse conduction band effective mass in silicon, and 
as = 5.3 X 10~^^ m is the Bohr radius for the hydrogen 
atom.^— We take the binding energy of a single electron 
to a single donor to be the experimentally observed value 
of -45.5 meV^a 

With these assumptions, the energy levels of the double 
donor, single electron system can be determined by the 
same variational linear combination of atomic orbitals 
(LCAO) technique used to calculate the eigenvalues of 
an moleculei^ In contrast to a real molecule, 
however, the position of the donors is fixed within the 
silicon lattice, and so it is not necessary to minimize the 
energy with respect to the donor separation. We plot 
the tunnelling frequency. A, as a function of the donor 
separation, i?, in Fig. [21 



III. OUTLINE OF THE SCHEME 

The scheme for measuring the decoherence rate for the 
charge qubit proceeds as follows. Initially, the electron 
is prepared in the \L) state by placing a large electric 
field across the double donor system, such that the bias 
term in the qubit Hamiltonian takes the value eg 3> A. 
£o must be chosen such that the total energy gap for 
the two level system satisfies HE = h^y A'^ + £q > ksT. 
The electron will then relax to the ground state, which 
is strongly localized on the left donor site. 

The system is then placed in the symmetric superpo- 
sition state, -I- |i?))/\/2, by adiabatically sweeping 
the bias field to the symmetry point s{t) = 0. The bias 
field sweep should be performed quickly, so that there 
is negligible decoherence during the sweep. However the 
sweep must not be made too quickly, or there will be 
coherent non-adiabatic transitions into the excited state. 
We discuss these requirements in more detail in Sec. 

The bias field is held at zero for a time ihoW- During 
this time, as a result of the interaction with the environ- 
ment, the qubit will lose coherence. This loss of coherence 
will be manifest in the decay of the off diagonal elements 
of the qubit density matrix (written in the \L), \R) basis) 
at a rate Fq. 

The bias field is then swept, rapidly but adiabatically, 
to a large negative value, — eo, and held at this value while 
the position of the electron is read out by a nearby elec- 
trometer. Thus e{t) has the following time dependence: 
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e(t) = tsw < i < isw + ^hold 

^hold 

(3) 

where tgw is the time taken for each bias sweep. 

Provided the electrometer can determine the position 
of the electron charge in a timescale shorter than the re- 
laxation time for the qubit under the large bias (— eo)i 
the readout process will correspond to a strong quantum 
measurement in the {\L), \R)} basis. As we discuss fur- 
ther in Sec. 0this measurement can be implemented by 
existing radio frequency single electron transistor (RF- 
SET) technology, 

By repeating the above preparation, bias sweep, and 
measurement steps a number of times, one can deter- 
mine the probability. Pa, of finding the electron on the 
right donor site at the end of the sweep. If the decay 
of the off-diagonal elements of the density matrix (in the 
\L), \R) basis) during the time thoid for which the bias is 
held at e = is negligible, then the electron will coher- 
ently tunnel through to the right donor site as the bias 
is swept through to e = — eo- The final state will be, 
approximately, the pure state, \R). Thus the observed 
probability of finding the electron on the right donor will 
be close to unity. Conversely, if there has been substan- 
tial decay of the off-diagonal elements during the 'hold' 
part of the evolution, the final state will be mixed, and 
the observed Pr will be substantially less than unity. 
The off-diagonal density matrix elements are expected to 
decay over a timescale ^q^- Thus repeating the whole 
procedure for different values of thoid , and plotting Pr as 
a function of ihoid should allow one to determine Fr^. 



IV. ESTIMATING THE DECOHERENCE RATE 
BY ADIABATIC TUNNELLING 

A. The model of decoherence 

In order to study the effects of the environment on the 
qubit, we model the environment as a bath of harmonic 
oscillator modes, linearly coupled to the z— component of 
the qubit, via the spin boson Hamiltonian, 



(4) 



-ffqb is the qubit Hamiltonian of Eq. . The second term 
describes the coupling between the position degree of 
freedom of the electron (cr^), and the displacement op- 
erators for the bath modes (a| + Qi), where the A^'s are 
coupling coefficients. The last term represents the free 
Hamiltonian of the oscillator bath, where the oj^'s are 
the angular frequencies of the bath modes. 



The spin boson Hamiltonian has been studied 
extensivelyji^ii^ The behavior of the system depends cru- 
cially on the spectral density of the bath, defined as 



J(^) 



XfS{uj - LOi) . 



(5) 



In general, the dynamics of the spin-boson model, for an 
arbitrary spectral density, is rather complicated. For the 
purposes of this work, however, a number of simplifying 
assumptions can be made. Firstly, we assume that the 
sweep of the bias field, e{t), is made sufficiently slowly 
that an adiabatic approximation can be employed. In 
particular, we require^^ 



(We discuss this adiabaticity requirement further in Sec. 
IVl ) Secondly, we assume a weak system-bath coupling, 
such that J{k^T/h) <C E{t), throughout the sweep. Fi- 
nally, we take the initial state of the qubit to be the 
thermal state 

_ exp(-gqb.o/fcar) , , 

^° " tr[exp(~i/qb,oABT)] ' ^'^ 

where -ffqb.o is the initial qubit Hamiltonian, i.e. Eqn. 
Q with e{t) — £o- Note that po is diagonal in the en- 
ergy eigenbasis of the initial qubit Hamiltonian. Under 
these assumptions, the density matrix of the qubit is al- 
ways diagonal in the instantaneous energy eigenbasis of 
the qubit Hamiltoniani2& In this case, the Bloch vector 
^(^) — ii'^x), (cy), (cz)) always lies parallel to the vector 
B — (A, 0, e{t)), and the dynamics can be understood by 
considering the evolution of r{t) — \r{t)\, the length of 
the Bloch vector. The evolution of r{t) under the above 
assumptions, is given by^ 



f(<) = -r(t) (rit) - req(O) 



(8) 



where the instantaneous relaxation rate T(t) depends on 
the spin boson model parameters ji^** 



nt) = ^smH J {E{t)) coth(^^^ 



(9) 



where 9 — tan^^(A/£). rcq{t) is the thermal equilibrium 
value of the Bloch vector, evaluated for the instantaneous 
energy gap of the system, r^qit) = tanh{hE{t)/2kBT). 

At low frequencies, the spectral density of the bath typ- 
ically has a power law behavior fi^sii J(w) cx w*, where 
the exponent s depends on the nature of the environment. 
Two potentially serious sources of decoherence in this 
system are a deformation potential coupling between the 
qubit and acoustic phonons, and an electrostatic coupling 
to Nyquist-Johnson voltage fluctuations, which may orig- 
inate in the surface electrodes used to control the qubit 
Hamiltonian parameters. The former is described by a 
superohmic spectral density (s > 1). However, as we 
show in Appendix ^ with a judicious choice of donor 
configuration, the decoherence rate due to phonons can 
be made negligibly small, and therefore we neglect it in 
what follows. 
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B. Results for Ohmic damping 

In this section, we concentrate on the case of deco- 
herence due to Nyquist-Johnson voltage noise, which is 
characterized by a bath with an Ohmic spectral density 
(s = 1). At low frequencies, the spectral density may be 
writter 



,13.14 



J{uj) — 2auj , 



(10) 



where a is a dimensionless parameter which characterizes 
the strength of the system-bath coupling. 

In order to estimate a, we follow a similar procedure 
to that employed in Ref. |3 We first define the bath 
operator 



X = 



(11) 
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FIG. 3: Sample evolution of the Bloch vector components 
and Tz, for part of the bias sweep. Also shown is the time 
profile of the bias sweep itself, e(t) (right axis). 



which couples to the Gz operator of the qubit, via the 
second term in Eq. Q . To proceed, we calculate the spec- 
trum of fluctuations in X in terms of the spectral density 
J{uj), and relate this to the spectrum of Nyquist-Johnson 
fluctuations in the surface gates. For a bath of harmonic 
oscillator modes in thermal equilibrium at temperature 
T, the Fourier transform of the symmetrized correlation 
function of this operator takes the form 



Sxiuj) 



i([X(t + r),X(t)]^)e— dr 



= ttJ{uj) coth 



2kBT 



(12) 



where [A, B]^ = AB + BA denotes an anti-commutator, 
X{t) = e-^Ht/hj^^-iHt/h ig ^Yie bath operator in the 
Heisenberg picture, and (O) = tr[Opcnv] denotes the ex- 
pectation of O for an environment in a thermal equilib- 
rium state, Pcnv 

For noise due to voltage fluctuations, X may be re- 
lated to a perturbation SVlr in the potential difference 
between the two donor sites by 



X 



eSVi 



LR 



2h 



(13) 



SVlb^ is related to the voltage fluctuations in the surface 
gates by 



SV, 



LR 



(14) 



where the dimensionless parameter /3 quantifies the elec- 
trostatic coupling between the surface gates and the 
donor sites, and is determined by the device geometry. 
For the geometry shown in Fig. ^ /3 may be approxi- 
mated by elementary electrostatics as 



2 In 



( 12. 



In 



d-rg 



In 



(d 



(15) 



where d is the distance between the two surface elec- 
trodes, a is the thickness of the oxide layer, ri is the 
distance between the left donor and the left electrode, 
r2 is the distance between the left donor and the right 
electrode, ro is the effective radius of the electrode, and 
ei and £2 are the dielectric constants of the oxide and 
silicon layers, respectively. In deriving this expression 
we have assumed that the gates may be represented by 
long, cylindrical conductors, and that tq <ti a, d. Using 
the values for ri , r2 , d, and a given in Fig. ^ and taking 
rQ — 5 nm, ei = 4, and £2 — 12, we find f3 = 0.17. 

Substituting Ea. (|14|l into Ea. ljl^ll . and calculating the 
corresponding power spectrum yields 



S'x(t^) 



4ft2 



Svito) ■ 



(16) 



For Nyquist-Johnson noise, the voltage fiuctuations are 
characterized by^^ 



Sviuj) 



= ^gatefi^COth(^) 



(17) 



where i?gatc is the impedance of the circuit which gener- 
ates the gate voltages, and T is the corresponding noise 
temperature. Substituting this expression into Eq. H16() 
and comparing with Eqs. (|10|l and l|12|) . we find that the 
system-bath coupling parameter is 



gate 



(18) 



where Rq — hje^ — 25.8 kJl is the quantum resistance. 
Taking i?gatc = 50ri and (3 ~ 0.17, we have a ~ 1.4 x 
10-^ 

We numerically solved Eqs. ijHJl and lO for the Ohmic 
spectral density of Eq. H10|) . assuming a bias sweep e(t) 
of the form described in section ITTll Figure |3| shows the 
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FIG. 4: Probability of finding electron on the right donor, 
Pa, at the end of the bias sweep, as a function of fhoid- The 
broken line represents the reciprocal of the decoherence rate 
for zero bias, Tq^. 



evolution of the x and z components of the Bloch vector, 
r{t), for a bias sweep with parameters Eq = 5 x 10^^ s~^, 
isw — 10"'' s, and ^hoid = 10~* s. We assume also that 
A ^ 10i°s-i and T = lOK. 

Figure^ shows the resultant probability, Pr = (1 — 
rz)/'^, that the electron is found on the right donor at 
the end of the sweep, as a function of thoid- The other 
parameters used in this calculation are the same as those 
used in Fig. 13 For values of ^hoW ^ Fq"^, Pr is close 
to unity, indicating that the electron has coherently tun- 
nelled from the left donor site to the right donor site. 
Note that Pr saturates to a value slightly less than unity, 
as a result of a small amount of decoherence during the 
'sweep' parts of the evolution. For values of thoid ^ Tq 
the resultant probability is substantially less than unity, 
indicating a loss of coherence during the 'hold' part of 
the evolution, due to interaction with the environment. 
The transition between these regimes occurs at a value 



of thoid ~ Fg . Thus measuring Pr at the end of the 
sweep provides a method for estimating the decoherence 
time Fo, and hence for estimating the strength of the 
system-environment coupling. 



V. DISCUSSION 

In order that the transition from coherent tunnelling 
[Pr « 1 at the end of the sweep) to incoherent behavior 
[Pr substantially less than 1) can be observed, and hence 
Fq be determined, the parameters A, Eq and tgw must 
satisfy a number of constraints. Firstly, we require that 
at time t = Q, the electron must be strongly localized 
at the left hand donor site. This can be achieved by 
placing a large bias Eq across the double donor system, 
and waiting for the donor to relax to its ground state. 
This implies that we require 



Secondly, we require that the minimum energy gap be- 
tween the ground and excited states satisfies 



(20) 



otherwise the system will simply remain in its ground 
state, throughout the bias sweep, and it will not be pos- 
sible to observe the effects of decoherence. 

Thirdly, coherent, non-adiabatic transitions into the 
excited level should be minimized. The problem of non- 
adiabatic transitions in two level systems was considered 
by Landau^^ and ZenerSS,. The results of Ref. '25* are 
directly applicable to the present work. For negligible 
nonadiabatic transitions, we require 



£0 



> 1. 



(21) 



Finally, it is necessary to ensure that relaxation at the 
end of the bias sweep (when e{t) < —kBT) is negligible. 
If there is significant relaxation over the last part of the 
bias sweep or during the measurement process, the sys- 
tem will be found to be in its ground state regardless of 
ihoid, and it will not be possible to observe the effects of 
decoherence. The probability that the electron will relax 
into the ground state, over the last part of the sweep, is 
approximately 



elax 



T{t)dt . 



(22) 



where F(t)is the relaxation rate of Eq. and is 

the time for which e{t) = e* = —ksT, and t/ — 2isw + 
thoid is the time corresponding to the end of the sweep. 
Performing the integral, and requiring that Prciax ^ 1, 
we have 



aTrA^ 

£0 



ln(^ 



< t: 



(23) 



In arriving at this expression, we have made the approx- 
imation coth{hE{t)/2kBT) Ri 1 for t > tt,, and that 
£o > E*. 

In order that there is no significant relaxation to the 
ground state during the measurement process, we require 



£0 



(24) 



where tmcas is the characteristic time for the electrometer 
to detect the presence or absence of the electron on the 
right donor site. We take tmcas — ^fJ-s, which is readily 
achievable with existing radio frequency single electron 
transistor (RF-SET) technology»i2.Sfl*Sl 

In the preceding analysis, T corresponds to the noise 
temperature of the electronics which generate the bias 
sweep. Taking T = 10 K, we find that kBT/h = 1.3 x 
10^^ s~^. In order to satisfy the inequalities of Eq. I|19|) 

10 , 



(19) and Eq. we choose A = 10^" s'^ and Eq = 5 
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10"'^^ s~^. The inequality of Eq. H21|l can then be satisfied 
if we choose tgw — 10^^ s. With these parameters, Eq. 
(I23|l and Eq. 124(1 imply that unwanted relaxation is 
negligible provided a < 3 x 10^'^. Comparison with our 
earlier estimate, from Sec. IIVBI of a « 10^^, suggests 
that the experiment is indeed feasible. 

A central element of the scheme introduced in Sec. IIIII 
is that the bias field, e(t), is held at zero for a time tjjoid- 
This bias field will be related to V^ato , the voltage across 
the surface electrodes in Fig. However, imperfections 
in the fabrication of a real device, and the existence of 
other surface electrodes (for instance, the plunger gate 
used to tune the RF-SET), may alter the potential land- 
scape in the vicinity of the donors, leading to a small 
systematic error, Se, in the bias field. This will lead to 
a systematic error in the observed value of the decoher- 
ence rate. According to Eq. (O, provided Me ^ fcsT, 
the observed rate will be 



where Fq is the decoherence rate evaluated for & = 0. 
Thus the true decoherence rate can be inferred by de- 
termining Fg for a range of different offset voltages and 
fitting the results to Eq. H25|l . Note that for sufficiently 
small offsets, 5e < A, we have 

To « (l - |J) To , (26) 

i.e. the error in the observed decoherence rate is only 
quadratic in the offset error. 

In our discussion of decoherence mechanisms, we have 
not explicitly considered errors due to background charge 
fiuctuations. These fluctuations vary from sample to 
sample, and typically have a 1// spectrum with a shoul- 
der at 100 — 1000 Hzi2P This timescale is longer than the 
time taken for each preparation, sweep, and measure cy- 
cle. Background charge fluctuations will therefore have 
the same effect as adding a small random offset bias, 
de, which may vary between cycles, but will be essen- 
tially constant over each bias sweep. As described above, 
the effect of such an offset becomes unimportant pro- 
vided that, in a given sample, the charge fluctuations are 
sufficiently small that the corresponding offsets satisfy 
Se < A. 



VI. CONCLUSION 

In summary, we have proposed and analyzed, theo- 
retically, an experimentally feasible scheme for directly 
determining the decoherence rate for a solid state charge 
qubit consisting of a single electron bound to a pair of 
donor impurities in a semiconductor host. The qubit 
is manipulated by adiabatically varying the bias term 
in the Hamiltonian. For a specific implementation us- 
ing phosphorous donors embedded in a silicon host, we 



have theoretically obtained quantitative estimates for the 
Hamiltonian parameters, and for decoherence rates cor- 
responding to interactions with both acoustic phonons 
and voltage fluctuations. We have analyzed various con- 
straints which must be satisfied in order for the exper- 
iment to be feasible. We have found appropriate, ex- 
perimentally achievable parameters which satisfy these 
constraints. Our results indicate that the control field 
needs to be manipulated with a time resolution of tens 
of nanoseconds, which is well within reach of current 
technology. Performing this experiment would be a vital 
step towards the implementation of a scalable solid state 
quantum computer. 
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APPENDIX A: RELAXATION DUE TO 
PHONONS 

In this appendix, we estimate the decoherence rate 
due to interaction with acoustic phonons. The prob- 
lem of electron scattering by acoustic phonons in sili- 
con was originally considered by Bardeen and Shockley^ 
More recently, electron relaxation, due to phonons, 
in low dimensional semiconductor systems has been 
considered f^^^ Due to the confinement of the electrons 
in these systems, and the resulting discrete spectrum of 
the electronic energy levels, relaxation due to phonons is 
suppressed. 

The rate for phonon emission in confined systems 

^31. 32 

X J dn,\{i^f\e^^"\i;,)\\ (Al) 

where D is a deformation potential, p is the density 
of silicon, Cg is the speed of sound, HE is the energy 
difference between the initial and final electron states, 
nB(w,Tph) = (exp(?ia;/fcB7ph) — 1)~^ is the Bose occu- 
pation function for a bath of phonons at temperature Tph, 
and Qi f is the wave number of the emitted phonon. qi f is 
fixed by the energy gap between the ground and excited 
states and the phonon dispersion relation as qif = E/cs- 
The integral in Eq. I|A1|I is over all solid angles in q 
space, and is evaluated for q — qif. In general, owing 
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FIG. 5: Approximate relaxation rate due to LA phonons, 
Fph, as a function of donor separation, for zero bias (e = 0). 
We assume Tph = 0.1 K, D = 3.3 eV, p = 2.33 kg m"^ 
Cs = 9.0 X lO'^m s ~^ and ao — 3 nm. 



to the anisotropy of the crystal, both D and Cg wih be 
tensors. However, for the purpose of gaining an order of 
magnitude estimate of Fph, we will ignore these subtleties 
and treat these quantities as being isotropic. 
The initial and final electron states are 



= cos-|L) +sin-|i?) 



1^/) =sin-|i) -cos-|i?) 



(A2) 



(A3) 



where 6* = tan ^{A/e), and \L) = {nag) 2 exp(— ra/as) 
and \R) = (7ra^)~2 exp(— rb/as) are Is-orbitals, with 
Bohr radius Ub, localized on the left and right donor site 
respectively. Taking the origin to be the mid point of the 
line joining the two donors, we have = r + ■^Rux and 
fb = r — ^Rux, where Ux is the unit vector along the line 
joining the two donors, and R is the donor separation. In 
these coordinates, the matrix element in Eq. (|Aip may 
be written 



sma 



dVf 



-2rb/aB 



+0((i?/as)'e-^/"") 



(A4) 



The last term in this expression may be neglected for 
donor separations R 3> 3as. Performing this integral, 
with the aid of the convolution theorem, we find that 



-16i sin sin((7j,i?/2) 



(A5) 
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FIG. 6: Approximate relaxation rate due to LA phonons, 
Fph, as a function of bias energy, e, for donor separations 
R = 35.0 nm, R — 37.5 nm, and R = 40.0 nm. We assume 
Tph = 0.1 K, D = 3.3 eV, p = 2.33 kg m-^ = 9.0 x 
lO^m s and ao — 3 nm. 



where Qx is the component of the phonon wavevector 
along the line joining the two donors. Substituting this 
expression into Eq. (|A1|I . and performing the integral 
over all solid angles, gives 



1] [l-sinc(g,/i?)] 



nphcs [(qifaB)^ + 4] 



(A6) 

Note that this rate is, in general, a function of the lattice 
temperature, Tph, the distance between the donors, R 
(which fixes A as shown in Fig. and the bias between 
the donors, s. 

In Fig. |31 we plot Fph as a function of donor separa- 
tion for zero bias (e = 0). In Fig. |Blwe plot Fph for a 
nonzero bias, for three different donor separations (i? = 
35.0 nm, R = 37.5 nm, and R = 40.0 nm). In these cal- 
culations we assume Tph = 0.1 K, p — 2.33 kg m~'^, 
D = 3.3 eV, Cs = 9.0 x 10^ m s'^ (LA phonons) and 
flo = 3 nmi2^ Note that the lattice temperature Tph used 
here is much less than the effective noise temperature 
(T — 10 K) assumed in Sec. the latter is due to noise 
in the electronics used to generate the bias sweep, which 
is typically much greater than the sample base tempera- 
ture. 

These results indicate that Fph is a strongly decreas- 
ing function of the donor separation (for i? > 20 nm). 
For donor separations greater than about 35 nm, Fpj^^ is 
significantly longer than the other relevant timescales in 
the problem, and therefore we are justified in neglecting 
phonons as a source of decoherence in this system. 
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